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ABSTRACT 

We analyzed RXTE-PCA observations of a recent outburst of the X-ray pulsar 
XMMU J054134. 7-682550. We calculated the pulse frequency history of the source. 
We found no sign of a binary companion. The source spins up when the X-ray flux 
is higher, with a correlation between the spin-up rate and X-ray flux, which may be 
interpreted as a sign of an accretion disk. On the other hand, the source was found 
to have an almost constant spin frequency when the X-ray flux is lower without any 
clear sign of a spin-down episode. The decrease in pulsed fraction with decreasing 
X-ray flux was intrepreted as a sign of accretion geometry change, but we did not find 
any evidence of a transition from accretor to propeller regimes. The source was found 
to have variable pulse profiles. Two peaks in pulse profiles were usually observed. We 
studied the X-ray spectral evolution of the source throughout the observation. Pulse 
phase resolved analysis does not provide any further evidence for a cyclotron line, but 
may suggest a slight variation of intensity and width of the 6.4 keV iron line with 
phase. 

Key words: X-rays: binaries, pulsars: individuahXMMU J054134. 7— 682550 , stars: 
neutron, accretion, accretion discs 



1 INTRODUCTION 

The X-ray pulsar XMMU J054134.7-682550 in the Large 
Magellanic Cloud (LMC) was discovered as a High Mass 
X-ray Binary (HMXB) in an XMM-Newton survey of the 
LMC (Shtykovskiy & Gilfanov 2005). The source was then 
listed as an HMXB in a catalog of HMXBs in the Magellanic 
Clouds (Liu, van Paradijs & van den Heuvel 2005). 

From a routine scanning of Swift-BAT data on 3 August 
2007, XMMU J054134.7-682550 was found to be in out- 
burst with an increase in flux levels from ~ 10 — 20 mCrab 
to ~ 50 mCrab (Palmer, Grupe & Krimm 2007). Spectral 
features, the transient nature, and pulsations at 61.6 s de- 
tected in RXTE observations during the outburst were con- 
sidered to be an indication of the high mass X-ray binary 
nature of the source (Markwardt, Swank & Corbet 2007). 
These authors also noted cyclotron absorption features at 
~ 10 and 20 keV in the X-ray spectrum of the source. 

In this paper we investigate the X-ray timing 
and spectral properties of the recent outburst of 
XMMU J054134.7-682550 from 9 August 2007 through 
1 November 2007 (MJD 54321 - MJD 54405). During 
this period the 3-20 keV flux declined from ^ 3.3 x 
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10"^"ergs"^cm"^ to ~ 6 x 10^^^ ergs"^cm" . We track 
the changes in spin period of the source, relating these to 
changes in the accretion geometry, and search for signatures 
of binarity in the lightcurve. We also study the evolution 
and flux dependence of the pulsed fraction through the out- 
burst. Furthermore, we explore the spectral evolution of the 
source over the outburst and with neutron star spin phase 
and investigate the possibility of cyclotron features in the 
spectrum. 



2 INSTRUMENTS AND OBSERVATIONS 

The Proportional Counter Array (PGA, Jahoda et al. 2006) 
on board the Rossi X-ray Timing Explorer (RXTE) con- 
sists of five identical proportional counters coaligned to the 
same point in the sky, with a total geometric area of ap- 
proximately 6250 cm^, a field of view to 1 square degree, 
and operating in the 2-60 keV energy range. In May 2000 
and in December 2006, PCUO and PCUl, respectively, lost 
their propane veto layers. Without a propane veto layer, 
these PCUs have higher background levels. We have ex- 
cluded the data from these PCUs in our spectral analysis. 
To obtain better energy channel counting statistics, we used 
only the data obtained from the top PCU layers for the 
spectral analysis. For both spectral and timing analysis, the 
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Table 1. PCA observations of XMMU J054134.7-682550. The 



MJD column refers to the start of the 


Obs ID 


Time 


Exposure 


93413-01- 


(MJD) 


(ks) 


01-00 


54321.108 


3.5 


03-00 


54327.014 


2.0 


04-00 


54330.669 


1.4 


04-01 


54330.733 


1.7 


05-00 


54333.796 


3.9 


06-00 


54336.954 


2.5 


07-00 


54339.707 


3.2 


08-01 


54341.875 


2.1 


08-00 


54342.197 


0.9 


09-00 


54345.526 


1.4 


09-01 


54345.798 


2.0 


10-00 


54348.613 


1.4 


10-01 


54348.678 


1.6 


11-00 


54351.595 


3.8 


12-00 


54354.950 


5.5 


13-00 


54357.048 


3.1 


14-00 


54360.133 


1.8 


15-00 


54363.926 


3.1 


16-00 


54366.932 


2.9 


17-00 


54369.753 


2.3 


18-00 


54372.491 


1.4 


18-01 


54372.560 


1.5 


19-00 


54375.628 


2.6 


20-00 


54378.914 


1.6 


20-01 


54379.121 


1.4 


21-00 


54381.995 


1.4 


21-01 


54382.06 


1.5 


22-00 


54384.42 


1.3 


22-01 


54384.474 


1.6 


23-00 


54387.420 


1.5 


23-01 


54387.486 


1.5 


24-00 


54390.619 


3.2 


25-01 


54393.519 


1.5 


25-00 


54393.712 


2.1 


26-00 


54396.704 


3.5 


27-00 


54399.067 


1.6 


27-01 


54399.143 


2.5 


28-00 


54402.451 


3.8 


29-00 


54405.807 


2.2 
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Figure 1. Long-term RXTE-ASM Lightcurve of 
XMMU J054134.7-682550. The inset shows a ~ 300 day 
portion around the outburst at ~ MJD 54365. Dots in the inset 
represent individual PCA observations listed in Table [l] 
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latest combined background models were used together with 
HEASOFT 6.3 to estimate the appropriate background. 

The long-term RXTE-ASM lightcurve of the 
source is presented in Fig. [1] RXTE observed 
XMMU J054134.7-682550 through a Type II outburst and 
28 observations were undertaken between 9 August 2007 
and 1 November 2007 with a total exposure of 127 ks (see 
Table A 3-25 keV background subtracted RXTE PCA 
lightcurve of the source is presented in Fig. [51 Although 
the number of active PCUs varied between one and three 
during the observations, count rates in Fig. [2] have been 
adjusted as if 5 PCUs were active using the correctlc tool 
in HEASOFT. 



3 TIMING ANALYSIS 

We used background corrected lightcurves for timing anal- 
ysis. These lightcurves were also corrected to the barycen- 



Figure 2. 61s binned 3-25 keV lightcurve of 
XMMU J054134.7-682550 



tre of the Solar system. From the periodogram around 61 s, 
a period was found giving the maximum x^- ^ template 
pulse profile from each RXTE observation was constructed 
by folding the data on this period. The template pulse pro- 
file comprising 20 phase bins was analytically represented 
by its Fourier harmonics (Deeter & Boynton 1985) and 
cross-correlated with the harmonic representation of aver- 
age pulse profiles obtained from each observation. The max- 
imum value of the cross-correlation is well-defined and does 
not depend on the phase binning of the pulses. We used a 
10-term unweighted harmonic series to cross-correlate the 
template pulse profile with the pulse profiles obtained from 
~ 400 — 500 s long segments of the observation. For each 
single or two succeeding RXTE observations, we obtained 
~ 5 — 12 arrival times. From the slopes of arrival times 
(5<j!> — 00 + Su{t — to), where 5(j> is the pulse phase offset 
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Table 2. Timing Solution of XMMU J054134.7-682550 



Parameter 
Epoch 

Spin Frequency(!/) 

Spin Frequency Derivative (i>) 

£> 

v 

V 

RMS Residual (pulse phase) 



Value 



MJD 54364.99681(4) 
1.6455121(7) X 10-2 Hz 
3.1412(1) X 10"" Hz s-i 
-1.2274(7) X 10-^'^ Hz s'^ 
-1.25(8) X 10-25 Hz s-3 
-8.0(1) X 10-31 Hz S-* 
0.056 



deduced from the pulse timing analysis, to is the mid-time 
of the observation, (j)o is the phase offset at to and Sv is the 
deviation from the mean pulse frequency), we corrected the 
pulse frequencies which are presented in Fig. [3] (top panel). 

We were able to phase connect all observation segments 
by using a fifth-degree polynomial in a Taylor expansion, 

5(t> =00 + Su{t - to) + ^i>it - tof + ^i>{t - tof 

L-O'U^tot + — vit-^tof (1) 
24 ^ ^120 

where v, i>, ii , v are the first, second, third and fourth 
derivatives of the pulse frequency of the source. Pulse resid- 
uals from the subraction of the fifth degree polynomial are 
presented in Fig. [S] 

Table [2] presents the timing solution of 
XMMU J054134.7-682550. In the top panel of Fig. H 
we present independent pulse frequency measurements 
for each observation and best fit model given by timing 
solution. In the middle and bottom panels, pulse frequency 
and pulse phase residuals are presented. In both residual 
panels, we have not seen any signature of a Doppler shift 
due to binary orbital modulation. 

Folding 3-8 keV, 8-13 keV, 13-18 keV and 18-25 keV 
light-curves of each observation with the pulse frequencies 
presented in Fig. [31 we obtained energy dependent pulse 
profiles of the source. We present 12 sample pulse pro- 
files from three observations in Fig. U We also obtained 
pulse profiles from 3-25 keV light-curves. Using these pulse 
profiles, we calculated pulsed fractions using the definition 
(-Fmax - -Fmin)/(-Fmax + -Fmin)) where F,nax and F,nin are the 
highest and lowest fluxes of the phase bins. In Fig. [S] we 
present evolution and X-ray flux dependence of the pulsed 
fraction of the source. X-ray flux values are obtained using 
the spectral analysis of the source (see next section) . 

For each sequence of 3-4 pulse frequency measurements, 
we fit a straight line to the pulse frequency time series in 
order to find the average pulse frequency derivatives. We 
plot the frequency derivative as a function of the average 
3-20 keV X-ray fiux in Fig. EJ 



4 SPECTRAL ANALYSIS 

Spectrum, background and response files for the observa- 
tions were created using HEASOFT 6.3. We did not use 
channels corresponding to photon energies below 3 keV due 
to uncertainties in the background at these energies; neither 
did we use channels above 20keV (25keV for the brightest 
observation on MJD 54321) due to low count statistics. No 



T 



\ 

i 



J I I L 

10 

.3-20 keV Unabsorbed Flux (10-"erg/s} 

Figure 6. Variation of frequency derivative with X-ray fiux. 

systematic error was included in the fits, as statistical errors 
dominate for these data. 

The brightest observation, on MJD 54321, was fit with 
an absorbed cut-off power law (White, Swank & Holt 1983) 
and a Gaussian emission component to account for a weak 
iron line at 6.4 keV. The neutral hydrogen density was un- 
constrained by the fit, and thus fixed to the line-of-sight 
value of 2.5 x lO^^ cm-^ (Kalberla et al. 2005) . The fit, 
along with Ax'^ residuals is shown in Fig. [7] and reported in 
column 2 of Table |3] 

Given that this model does not provide an optimal de- 
scription of the observed spectrum (xt ~ 1-7), and that 
Markwardt et al. (2007) reported possible cyclotron line fea- 
tures at 10 and 20 keV, we introduced cyclotron absorption 
features into the model at these energies. The model fit was 
clearly improved, however, the F-test probability of this im- 
provement being by chance is quite significant: 0.018. 

An alternative continuum model was then used ( 
powerlaw * highecut) which provided a better fit (see Ta- 
ble [3] and residuals in Fig. [7| of the continuum without the 
need for a cyclotron feature (see Fig.O, lower panel). In this 
fit we have allowed the neutral hydrogen density to vary, but 
it is still poorly constrained by the model. This model intro- 
duces a discontinuity in the spectrum at the cutoff energy, 
and such a discontinuity can disguise the presence of a cy- 
clotron feature. Thus, our spectral analysis cannot confirm 
or rule out a cyclotron line in this source. 

In order to investigate the variability of the spectrum 
through the outburst of XMMU J054134.7-682550, the best 
fit model above (power*highecut) was fit to all observations 
in Table [T] in 3-20 keV band using a fixed uh value of 8 x 
lO^" cm-2. For most observations an Fe line at ~ 6.4 keV 
was required to fit the spectra. From the spectral fits cut- 
off energy values were found to be between ~ 11 and ~ 21 
keV, and -Efoid energy values were found to be between ^ 4 
and - 33 keV. In Fig. [H we plot 3-20 keV unabsorbed X- 
ray flux, spectral parameters and reduced obtained from 
spectral fits as a function of observation days. 

Assuming a distance of 49kpc (Catelan, Cortes, 2008) 
and using the unabsorbed X-ray flux values plotted in Fig. [8] 
and listed in Table [TJ we found that X-ray luminosity of the 
source in the 3-20 keV band varies between ~ 2 x lO^'' - 9 x 
10^'' ergs-^ 
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Figure 4. Energy dependent pulse profiles for observations on MJD 54321.1 (top left panel), MJD 54351.6 (top right panel), and MJD 
54387.4 (bottom panel). The unit of y-axis of each plot is normalized count from the specified energy interval. 
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Figure 5. Evolution and X-ray flux dependence of pulsed fraction for XMMU J054134. 7-682550. 
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Figure 7. The upper panel shows the cutoff power law model 
fitted to the spectrum. The middle panel shows the Ax'^ residuals 
from this model fit. A dip at around 10 keV is clearly shown, but 
the model is a poor fit to the continuum in general. The lower 
panel shows the Ax'^ residuals from fitting a power law model 
with a high energy cutoff. The dip at 10 keV is no longer apparent, 
showing that a cyclotron component at 10 keV is not necessary. 



Table 3. Results of spectral analysis. Errors are 90% uncertain- 
ties. Flux is in ergcm" ^s^^ and the Fe normalisation in pho- 



Parameter 



cutoffpl power * highecut 



A^H X 10^1 cm-2 


2.5 (fixed) 


0.8lO;| 


r 


AO+O-Oi 
'^•^^-0.05 


0.95l0:0«9 


Scut (keV) 




14 7+0-8 


-Bfold (keV) 


14.310:5 




Efc (keV) 


fi 5+°-2 

D.O_Q 2 


O.O_0.4 


^Fc (keV) 


n .,-1-0.3 


1 0+0-5 


Fe normalisation 


2Atll X 10-4 


itl X 10-4 


Unabs Flux (3-25 keV) 


4.2 X 10-10 


4.2 X 10-10 


xl (d.o.f) 


1.7 (50) 


0.72 (48) 




?E 5 ^ i 5 $ . 



Time (days sirce MJD 54364,996811 

Figure 8. 3-20 keV unabsorbed flux in units of lO-^^ 
erg.s-l cm— ^, power law index, iron energy in units of keV and 
reduced x^ as a function of observation time. Errors indicate la 
confldence level. 



Fig. |9] shows the dependence of power law index on the 
flux. The power law index shows an anticorrelation with 
the flux, i.e. the spectrum gets softer with decreasing X- 
ray flux. Spectral softening accompanied with X-ray flux 
decrease can be considered to be a result mass accretion rate 
variations (Meszaros et al. 1983, Harding et al. 1984), and is 
also observed in other accretion powered X-ray pulsars like 
SAX J2103.5+4545 (Baykal et al 2002) and 2S 1417-62 
(Inam et al. 2004). 



5 PULSE PHASE RESOLVED 
SPECTROSCOPY 

As cyclotron features are strongly dependent on viewing an- 
gle, phase resolved spectroscopy was undertaken to further 
study the presence of possible cyclotron features. A sec- 
ondary objective was to observe whether any other spectral 
features changed with the phase of the neutron star spin. 
Because of the distance to the source and the weakness of 
the possible cyclotron feature, we focus here on the bright- 
est of the observations (01-00 at MJD 54321.1) to ensure the 
best signal-to-noise possible. 

The neutron star pulse period was arbitrarily divided 
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Figure 9. Variation of power law index with X-ray flux. 




trum, however, the xt was improved from 1.3 to 1.1 by the 
addition of a cyclotron feature at lOkeV. The F-test prob- 
abihty of this improvement being due to chance is 0.007, 
suggesting a rather marginal detection. 

To be consistent with the analysis of the phase average 
spectra, we also fitted an absorbed power law model with a 
high energy cutoff. Again, this proved similar to the phase 
averaged spectra in that the cutoff in the power law masked 
the residuals at lOkeV that were apparent in fitting the 
cutoff power law above. Once again, the hydrogen column 
density was fixed, and the energy of the iron line was fixed 
at 6.4 keV. With this model, too, the spectral slope showed 
little variation with pulse phase and no cyclotron features 
were required. Table |4] gives the parameters for the fits to 
the four spectra. 

The lower panel of Figure [TU] shows the variation in the 
width of the iron line for the four phases of our pulse pro- 
file for the high energy cutoff model with y-error bars the 
90% uncertainties from XSPEC. It is also clear from Table|3] 
that the iron line intensity follows the intensity of the source 
through the pulse profile, i.e. the iron line is strongest in the 
pulse peak. In Shoulder 1 the iron line is poorly constrained, 
and the low xl, together with the errors on the line width 
indicate that it is not strictly required in this phase. The 
errors are large due to the faintness of both the iron line 
and are probably infiuenced by uncertainty in the under- 
lying continuum. This makes it difficult to infer clearly any 
variation of the width of the line with phase from these data. 

Neither of the models fitted point strongly towards the 
presence of a cyclotron feature in this source, however, given 
the quality of the data, a cyclotron line cannot be ruled out. 
The data suggest that there may be some variation in both 
the width and intensity of the 6.4 keV iron line through neu- 
tron star spin phase, but further observations below 3keV, 
which will constrain the hydrogen column density, are re- 
quired to show any such variation with certainty. 



Figure 10. Upper Panel: Pulse profile for observation 01-00. The 
four phases studied are marked. Lower Panel: Width of iron line 
as a function of pulse phase using the power law model with a 
high energy cutoff. The y-errors are the 90% uncertainties given 
by XSPEC. 



into four phases according to the general morphology of the 
pulse profile, and in such as way that the chosen phases had 
broadly adequate signal-to-noise. The upper panel of Fig. 1101 
shows the pulse profile for observation 01-00 folded on a 
61.6 s period with the phase bins comprising two shoulder 
features and the rising and falling limbs of the main peak. 

Using the spectral fits to the phase averaged spectrum 
in the previous section as a guide, we fitted a cutoff power 
law model to the phase resolved spectra to see if there were 
any significant deviations from the phase averaged contin- 
uum. The hydrogen column density was fixed to the stan- 
dard LMC value of 2.5 x lO^^cm'^ (Kalberla et al. 2005). 
A simple absorbed cutoff power law was an adequate fit 
{xt < 1.1) to the data for the Shoulder 1 and Main 1 phases, 
but an iron line at 6.4 keV (fixed) was required in the Main 
2 and Shoulder 2 spectra. There was found to be little vari- 
ation in the slope of the spectra, and all four showed weak 
signs of an absorption feature at lOkeV, much the same as 
the phase averaged spectrum. In the case of the Main 2 spec- 



6 DISCUSSION AND CONCLUSION 

The lack of orbital signature in pulse frequency and phase 
residuals is evident in Fig. [3] We also searched for a possible 
periodicity in the RXTE-ASM lightcurve as a sign of the bi- 
nary orbit, but we did not find a sign of orbital modulation. 
The lack of orbital signature may be due to the fact that the 
orbital period of the binary sytem is larger than the total 
observation time span of 74 days. Future X-ray monitor- 
ing observations are needed to detect the signs of the orbit 
of this system. 

The transient nature of the source may be a sign of a 
Be type companion, i.e. the enhancement in the accretion 
rate may be due to the presence of a temporary accretion 
disk that forms when the neutron star enters denser equa- 
torial stellar wind of the Be companion. Thus, the outburst 
is probably a type II outburst from a Be star system. From 
Fig. [1] this outburst is seen to be the brightest outburst 
from this source since the launch of RXTE. 

In Fig. [S] a spin-up of trend of the source is evident 
especially for the first ~ 40 days of the observations. In 
Fig. |S] it is seen that flux is correlated with spin-up rate 
of the source. From Fig, [S] it is evident that observations 
corresponding to the lowest flux values have spin rates con- 
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Table 4. Parameters from the powerlaw model with a high energy cutoff fitted to the four phases indicated in Fig lf 01 The centroid of 
the iron line was fixed at 6.4 keV and the Nh at 2.5 X fO'^^ cm~^ in all cases. Errors are 90% uncertainties. 



Parameter 


Shoulder f 


Main f 


Main 2 


Shoulder 2 


<TFe (keV) 


0.5 ±0.3 


f.0±0.3 


1.6 ± 0.2 


1.5 ±0.2 


/po {phcm~^s~'^) 


(4 ±2) X 10^4 


(9 ±4) X fO-4 


(1.5 ± 0.5) X 10"3 


(9.1 ±0.5) X 10-4 


r 


0.9f ±0.02 


0.86 ±0.03 


0.91 ± 0.01 


0.90 ±0.02 


xl (58 d.o.f) 


0.69 


0.8f 


0.92 


1.10 



sistent with zero, showing that source is either spinning-up 
or it is at an almost constant spin frequency phase with- 
out any clear evidence of a spin-down episode. If we as- 
sume that observed X-ray luminosity is proportional to the 
bolometric luminosity (or mass accretion rate), then spin- up 
rate (which is proportional to torque exerted on the neutron 
star) and X-ray flux correlation can be explained by accre- 
tion from accretion discs when the net torque is positive 
and of the order of the material torque (Ghosh & Lamb 
1979; Ghosh 1993). There are transient X-ray pulsar sys- 
tems which show correlation between spin-up rate and X-ray 
flux: EXO 2030+375 (Parmar, White & Stella 1989; Wilson 
et al. 2002), A 0535+26 (Finger, Wilson & Harmon 1996a; 
Bildsten et al. 1997), 2S 1845-024 (Finger et al. 1999), 
GRO J1744-28 (Bildsten et al. 1997), GRO J1750-27 
(Scott et al. 1997), XTE J1543+568 (in't Zand, Corbet & 
MarshaU 2001) and SAX J2103.5+4545 (Baykal, Stark & 
Swank 2002, Baykal et al. 2007), 2S 1417-62 (Finger et al. 
1996b, Inam et al. 2004). Thus, the positive correlation we 
show in Fig. |S] may be an indication of accretion via an ac- 
cretion disk. However, without any orbital correction, it is 
not possible to analyze the nature of this correlation quan- 
titatively, since the calculated pulse frequencies might not 
accurately correspond to the pulsar's spin frequencies. 

Decreasing pulsed fraction with decreasing X-ray flux 
(Fig. O and changing pulse morphology (Fig. |3| may in- 
dicate changes in accretion geometry and even be an in- 
dicator of a transition from disc to wind accretion, or an 
ongoing transition from accretor phase to propeller phase. 
This is similar to the decrease in pulsed fraction and pulse 
morphology changes in the low X-ray flux observations in 
2S 1417-62 (Inam et al. 2004). The evidence that observa- 
tions corresponding to the lowest flux values have spin rates 
consistent with zero (Fig. [S} may also be an indicator of 
accretion geometry changes in the source. 

The decline of the X-ray flux together with the change 
in pulse proflle shapes and pulsed fraction may be an indica- 
tion of a transition from accretor to propeller stage. When 
the propeUer stage sets in, the great majority of accreting 
matter cannot reach the neutron star. To test this hypothe- 
sis, we can compare the factor of decrease in the bolometric 
luminosity with the theoretical expectation for accretion to 
propeller stage, 

A = 170Mi{f J?g 'Pp/^ (2) 

where Mi. 4 is the mass of the neutron star in units of 
1.4Mq, Re is the radius of the neutron star in units of 10^ 
cm and Po is the spin period of the neutron star in units of 
seconds (Corbet 1996; Campana & Stella 2000; Campana et 
al. 2002). The factor A becomes ^ 2.6 x 10'' for a neutron 



star with mass ~ 1.4Mq, radius ~ 10 cm and spin period of 
~ 61 s. If we assume that the factor of decrease in bolometric 
luminosity is comparable to that of 3-20 keV luminosity, we 
find that our observed factor of luminosity change is ~ 50, 
which is much less than the theoretical expectation for a 
transition to propeller stage. It is not likely that we observe 
any transition to propeller stage. 

As seen in Fig. [J] pulse profiles show high variability. 
Except 3-8, 8-13 and 13-18 keV pulse profiles of the obser- 
vation on MJD 54321.1 and 18-25 keV pulse profile of the 
observation on MJD 54351.6, pulse profiles have two peaks 
seperated by ~ 0.5 phase. 

We found that both the phase average spectrum could 
be fit with either a cutoffpl model, with a possible weak 
cyclotron feature at 10 keV, or that it could be fit adequately 
using a power law with high energy cut-off model without 
the need for a cyclotron feature. We used the power law 
model with high energy cutoff to map the spectral evolution 
throughout all observations (see Fig. [S]). 

We also constructed pulse phase resolved spectra from 
the brightest observation. From this analysis (see Fig. llOf) . 
we did not find any conclusive evidence of a cyclotron ab- 
sorption feature, possibly due to low count statistics of the 
source. We found a very marignal variation in the iron line 
width and intensity which will need to be confirmed through 
further observations in the soft X-ray band. 
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